function y = Cauchy_rand(m,b)
u = rand();
y = m - (b ./ tan(pi .* u));
end